Rotationally-Invariant Exchange Interaction: The Case of Paramagnetic Iron 
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We present a generalization of the spin-fluctuation theory of magnetism which allows us to treat 
the full rotational invariance of the exchange interaction. The approach is formulated in terms 
of the local density approximation plus dynamical mean- field theory (LDA-(-DMFT), providing 
a systematic many-body treatment of the effect of spin-density fluctuations. This technique is 
employed to study the electronic and magnetic properties of paramagnetic a iron. Our result for 
the Curie temperature is in good agreement with experiment, while the calculations with the Ising- 
type exchange interaction yield almost twice overestimated value. 
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The theoretical description of metallic magnets, espe- 
cially those containing transition metals, is one of the 
central problems in condensed matter physics. Even to- 
day, in view of the great technological importance of such 
materials, a detailed understanding of their electronic, 
magnetic, and structural properties at finite tempera- 
tures remains problematic. This is mostly due to the 
presence of local magnetic moments above the magnetic 
ordering temperature which complicates the problem 
considerably and reduces the predictive power of first- 
principles calculations. Various properties of metallic 
magnets can be understood by using the spin-fluctuation 
theory [ij with its most general form based on a func- 
tional integral formulation This formulation was 
employed to describe the formation of local moments 
in paramagnetic metals 0] by reducing the many-body 
problem to a one-particle problem in a fluctuating ex- 
ternal magnetic field and then evaluating the functional 
average. Stimulated by these results, different analyt- 
ical and numerical methods have been developed, e^. 
the well-known quantum Monte Carlo techniques [j, |5| . 
By taking into account fluctuation corrections to the 
mean-field approximation, the spin-fluctuation theory 
have shown to provide a good qualitative description of 
the Curie- Weiss law behavior of magnetic susceptibility. 
However, applications of this technique to describe, e.g., 
the a-j phase transition in iron, do not lead to satis- 
factory results 0. In particular, it predicts the bcc-fcc 
phase transition to occur below the Curie temperature, 
Tc, while, in fact, this phase transition occurs 150 K 
above Tc- 

The LDA-hDMFT method 0, a combination of the 
ab initio local density approximation (LDA) of the den- 
sity functional theory and dynamical mean-field theory 
(DMFT), nowadays has become a state-of-the-art ap- 
proach for the calculation of the electronic and magnetic 
properties of correlated electron compounds Q. Appli- 
cations of the LDA-t-DMFT to study transition metal 



compounds have shown to give a good quantitative de- 
scription of localized as well as itinerant electron states 
These calculations predict the correct values of 
the local magnetic moment and magnetization, while the 
magnetic transition temperature turns out to be signifi- 
cantly overestimated. It was proposed that this is caused 
by the single-site nature of the DMFT approach, which 
is not able to capture the reduction of magnetic transi- 
tion temperature due to long wavelength spin waves [§] . 
Nevertheless, the substantial overestimation of the mag- 
netic transition temperature was also observed in other 
compounds, even in those where non-local spin fluctua- 
tions are negligible [l^. This implies that the approxi- 
mate form of the local Coulomb repulsion, restricted to 
the Ising-type exchange interaction, can be the reason of 
this systematic overestimation. The recently proposed 
continuous-time quantum Monte Carlo algorithms [l2j 
as well as some other quantum impurity solvers al- 
low one to treat the Coulomb interaction in its general 
form retaining rotational symmetry of spin. However, 
applications of these techniques so far have been limited 
to simple model systems due to the high computational 
costs. Therefore, only the density-density part of the 
Coulomb repulsion, that implies the Ising form of the ex- 
change interaction, has been employed in the most mate- 
rial specific calculations of correlated electron materials. 
However, as we will show below on the example of para- 
magnetic iron, the retaining of spin rotational symmetry 
is crucial for the correct description of the magnetic prop- 
erties. 

In this Letter, we present the spin-fluctuation theory 
of magnetism which is formulated in the framework of 
the LDA+DMFT method. The approach provides a sys- 
tematic treatment of the effect of local electronic corre- 
lations by reducing the many-body problem to the func- 
tional integral over a fluctuating magnetic field on an 
effective impurity. The spin-fluctuation theory is gener- 
alized by replacing a scalar fluctuating magnetic field to 
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a vector one. This allows one to take into account the 
full rotational invariance of the exchange interaction in- 
stead of the approximate Ising-type form. The proposed 
method is employed to study the electronic and magnetic 
properties of paramagnetic a iron, resulting in the Curie 
temperature value which is in good agreement with ex- 
periment. 

We start with the simple Hamiltonian of the Coulomb 
interaction in the following form 

i^Coul = ^ X! ^"^o-^liyCT + ^ X! ~ J)n^an„a, (1) 

where n^g. denotes the electron number operator with the 
spin (T (=t, i) at the orbital /.i. Using the total electron 
number operator, = "^^^ n^^^, and the z-projcction of 

the spin operator, Sz — J2fj,i^fj.t ^ the Hamilto- 

nian can be rewritten as 

i?Coui = ^UNiN - 1) + ^JN - JSl (2) 

where Lf — U ~ J/2 is the average value of the Coulomb 
interaction. This Hamiltonian represents the density- 
density part of the Coulomb interaction and contains 
the exchange interaction in the Ising-type form. To re- 
store the spin rotational symmetry, one should replace 
the z-projection of the spin operator, Sz, to the vector 

spin operator, S. Therefore, the Hamiltonian with the 
rotationally-invariant exchange interaction reads 

^Coui = ^UN{N - 1) + i JiV - jf. (3) 

Following the spin-fluctuation theory, we neglect the 
charge fluctuations but preserve the spin dynamics (mag- 
netic moment fluctuations) . Employing the static mean- 
field approximation for the first term in Eq. ([3]) and intro- 
ducing the double counting correction for the Coulomb 
interaction, the Hamiltonian of the system can be ex- 
pressed as 

H = i?LDA + Uina - ndo)N - jf, (4) 

where i?LDA is the LDA Hamiltonian, Ud = (N) is the 
average number of 3d electrons, and n^o is the LDA value 
for Hd- 

In the DMFT approach the lattice problem with 
the Hamiltonian ^ is mapped onto a quantum impu- 
rity model. Using the general form of the Hubbard- 
Stratonovich transformation [2[, the partition function 
can be expressed as a functional integral 

Z = l Di{r) exp[-^ J^' Pir)dr]ZiO, (5) 

where 

Z(0 =Tr{ T, exp[-/3ijLDA - l3U{nd - ndo)N 

+ 2cJ^'C{r)hr]}. (6) 



Here, T^ denotes the time ordering operator, /3 the 
inverse temperature, and c = ^J'KJjfi. Function ^(r) 
stands for an effective magnetic field resulting in the po- 
tential V(t) = 2c^(t)S. The functional integral over all 
fluctuating fields gives a solution of the impurity prob- 
lem. 

In the functional integral formulation of the conven- 
tional spin-fluctuation theory, the fluctuating magnetic 
field in Eq. ([5]) is considered to be scalar. The gener- 
alization to a vector field corresponds to the transition 
from Eq. ([2|) to Eq. (jS]) and allows one to take into ac- 
count the spin-rotational symmetry, thereby extending 
the theory from the Ising-type exchange interaction to 
the full rotationally-invariant one. Dividing the imag- 
inary time interval [0, /3] on L slices of length At and 
using the Trotter breakup for the exponential operator 
in Eq. (|6]), the partition function Z(^) for a given ^(t) 
can be expressed as 

Z(i) ^ Tr |t, J](cxp[-Ari/o] cxp[y(TO])| , (7) 

where i?o = -ffLDA + U^nd — nda)N is the ^-independent 
part of the Hamiltonian. These equations are similar 
to those of the Hirsch-Fye quantum Monte Carlo (HF- 
QMC) method Q. The partition function can be written 
as 

L 

z^Y. ^M-^Y. ^^'(^')] n dctic;! (0]. (8) 

However, instead of a single spin-fiip as in the HF-QMC 
method, here one should stochastically change the value 
of the field ^(t) for a random value of imaginary time. 
Due to the rotational symmetry of the exchange in- 
teraction, the interacting Green function becomes non- 
diagonal in spin indexes. The computational scheme 
where the partition function is calculated with an aux- 
iliary vector magnetic field is referred below as J-QMC. 
By taking into account only the z-component of the field, 
the approximate form of the local Coulomb interaction, 
limited to the Ising-type exchange interaction, is assumed 
(referred as J^-QMC). 

Elemental iron is one of the most famous itinerant- 
electron fcrroniagnets which exhibits localized moment 
behavior above the Curie temperature, Tc- Although 
various properties of the low-temperature ferromagnetic 
state of Fe can be understood within the density- 
functional theory applications of these techniques 
to describe the paramagnetic state do not lead to satis- 
factory results. Clearly, an overall understanding of the 
properties of iron requires a formalism which takes into 
account the existence of local magnetic moments above 
Tc [H. Recent apphcations of the LDA-f-DMFT have 
shown to provide a qualitatively correct description of the 
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FIG. 1. (Color online) Partial densities of states obtained by 
the J-QMC (top panel) and HF-QMC (bottom panel) calcu- 
lations within LDA+DMFT. The Fermi level is indicated by 
the vertical (gray) line at zero energy. Insets: imaginary parts 
of the self-energies. 



electronic, magnetic, and structural properties of para- 
magnetic iron However, a quantitative 



agreement has been achieved only in terms of the reduced 
temperature T/Tc, while the calculated Curie tempera- 
ture was found to be about twice larger than the ex- 
perimental value of 1043 K [l^. The average Coulomb 
interaction in the Fe 3d shell is considerably smaller than 
the bandwidth, showing no evidence for the formation of 
Hubbard bands in the spectral function. However, due to 
the strong exchange interaction [11|, the local magnetic 



moments are formed that is accompanied by the loss of 
coherence for the metaUic states due to the scattering 
of electrons on the fluctuating spins. The charge fluctua- 
tions in paramagnetic iron are of high frequency with the 
3d electrons being far from the localization limit. These 
arguments make iron an ideal candidate for our study. 

To calculate the electronic structure of paramagnetic 
a iron within the LDA, the tight-binding linear muffin- 
tin orbital (TB-LMTO) method was employed The 
low-energy Hamiltonian containing the 4s, 4p, and 3d 
states has been constructed with use of the TVth-order 
muffin-tin orbital (A^MTO) method [2l|. In our calcula- 



tions, we used the value of the screened Coulomb inter- 
action, U = 2.3 eV, and the value of Hund's exchange, 
J = 0.9 eV, which are consistent with the previous esti- 
mations [i,[ni[iil23i. 

In Fig. [T] we present the partial densities of states and 
the corresponding imaginary parts of the self-energies 
obtained by the LDA+DMFT at /3 = 10 eV~^ The 



FIG. 2. (Color online) Spin-spin correlation functions on 
the real and imaginary energy (inset) axes calculated by the 
HF-QMC and J-QMC within LDA-fDMFT. 



HF-QMC and J-QMC calculations give quantitatively 
similar results, reproducing the splitting in the density 
of states of the eg-orbitals near the Fermi level caused by 
exchange interaction [ll[ . The splitting in the density of 
states of the f2g-orbitals is found to occur in the LDA 
calculation and hence can be attributed to the band- 
structure effects. In both approaches, the self-energies 
for the orbitals remain Fermi-liquid-like, while the 
ones for the eg orbitals diverge at low frequencies. The 
latter indicates the formation of local magnetic moments 
that iniplies more incoherent eg states and itinerant t2g 
states [llj. We note that at low frequencies the self- 
energies obtained by the J-QMC method are close to 
those of the HF-QMC. This indicates that the physics 
near the Fermi level is dominated by the spin fluctua- 
tions while the charge fluctuations play a minor role. 

In Fig. [5] we show our results for the orbitally- resolved 
spin-spin correlation functions on the real and imaginary 
energy axes. In both approaches, the non-Fermi-liquid 
behavior of the eg electrons yields a pronounced peak at 
zero energy of the real energy axis indicating the presence 
of local magnetic moments. The satisfactory agreement 
of the results obtained by the HF-QMC and J-QMC 
methods suggests that the effect of the charge fluctua- 
tions, neglected within the J-QMC approach, is minor. 

To proceed further we compute the uniform magnetic 
susceptibility as a response to an external magnetic field. 
The temperature dependence of the inverse uniform mag- 
netic susceptibility obtained by the LDA-I-DMFT shows 
a linear behavior at high temperatures (Fig. [3]) . This in- 
dicates the presence of local magnetic moments and cor- 
responds to the Curie- Weiss law, = 3(T — Tq)/ ix\ff^ 
where Tc is the Curie temperature, jjLef / the effective lo- 
cal magnetic moment. The results of the least-square fit 
to the Curie- Weiss law are shown in Fig.[3]by the straight 
lines. It is clearly seen that the HF-QMC method, lim- 
ited to the Ising-type exchange interaction, overestimates 



4 



2 2 



1 1 1 
■ HF-QMC (r^ = 2050 K) 


li,ff=2^5- 


y -QMC (T^ = 1750 K) 




- ♦ 7-QMC (r^ = 850K) X 






^ieff=2-52/ 




^xp _ 
-C 


1043 K 

, 1 





TRR 80 and FOR 1346 is gratefully acknowledged. 



1000 2000 



4000 6000 
Temperature (K) 



8000 



FIG. 



3. (Color online) Temperature dependence of the 
inverse uniform magnetic susceptibility obtained by the 
LDA+DMFT. The straight lines depict the least-squares fit to 
the Curie- Weiss law. The experimental value of Tc = 1043 K 
is denoted by (black) arrow. The experimental value of the 
local magnetic moment is = 3.13 /^s [l9l |. 



the Curie temperature value almost twice. The Jz-QMC 
approach, which has the Ising-type exchange interac- 
tion, gives a slightly smaller value of the Tq than the 
HF-QMC. This confirms the validity of the static ap- 
proximation for the charge degrees of freedom. Taking 
into account the full rotationally-invariant exchange in- 
teraction, our calculations result in a substantial decrease 
of the Tc value, which is now found to be in satisfactory 
agreement with experiment. These findings are compati- 
ble with the results of the recent two-band model studies 

In conclusion, we presented a generalization of the spin 
fluctuation-theory of magnetism which allows one to take 
into account the full rotational invariance of the exchange 
interaction. The approach is formulated in terms of 
the LDA-I-DMFT method, providing a systematic many- 
body treatment of the effect of spin-density fluctuations. 
We employed this new technique to study the electronic 
and magnetic properties of a iron. Our results agree well 
with experiment and show that the overestimation of the 
Curie temperature by LDA-I-DMFT is mostly related to 
the approximate (Ising-type) treatment of the exchange 
Coulomb interaction rather than to the single-site nature 
of the DMFT. 
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